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Abstract 

We analyze the CP-violating electric dipole form factor of the nucleon in the framework of covariant baryon chiral 
perturbation theory. We give a new upper bound on the vacuum ang le, |0o| < 2.5 -10" 10 . The quark mass dependence 
of the electric dipole moment is discussed and compared to lattice QCD data. We also perform the matching between 
its representations in the three- and two-flavor theories. 
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1. The neutron electric dipole moment (nEDM) is a sensitive probe of CP violation in the Standard 
Model and beyond. The current experimental limit d n < 2.9 • 10~ 26 ecm [1] is still orders of magnitude 
larger than the Standard Model prediction due to weak interactions. However, in QCD the breaking of 
the U(1)a anomaly allows for strong CP violation, which is parameterized through the vacuum angle 9q. 
Therefore, an upper bound on d n allows to constrain the magnitude of 9q. New and on-going experiments 
with ultracold neutrons strive to improve these bounds even further, see e.g. [2] for a very recent review. On 
the theoretical side, first full lattice QCD calculations of the neutron electric dipole moment are becoming 
available [3-5] . These require a careful study of the quark mass dependence of the nEDM to connect to the 
physical light quark masses. In addition, CP-violating atomic effects can be sensitive to the nuclear Schiff 
moment, which receives a contribution from the radius of the nucleon electric dipole form factor, see e.g. [6]. 
It is thus of paramount interest to improve the existing calculations of these fundamental quantities in the 
framework of chiral perturbation theory. In [7] , the electric dipole moments of the neutron and the A were 
calculated within the framework of U(3)l X U(3)r heavy-baryon chiral perturbation theory and an estimate 
for #o was given (for earlier works utilizing chiral Lagrangians, see [8,9]). In [10], the electric dipole form 
factor of the nucleon was analyzed to leading one-loop accuracy in chiral SU(2), thus in that calculation 
the form factor originates entirely from the pion cloud. The strength of the form factor was shown to be 
proportional to a non-derivative, time- reversal- violating pion-nucleon coupling g^NN that could only be 
estimated from dimensional analysis. Furthermore, the leading contributions to the nEDM at finite volume 
and in partially-quenched calculations were considered in [11], and in [12] the leading order extrapolation 
formula using a mixed action chiral Lagrangian is given. In this Letter, we extend the results of [7,10] to 
higher order based on a covariant version of U(3)l X U(3)r baryon chiral perturbation theory. This allows 
to make contact to the lattice QCD results from [4] and by matching, we can also get more insights into the 
nucleon electric dipole form factor and the size of the coupling constant g-nNN ■ 
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2. The electromagnetic form factors of the nucleon (in the presence of P and CP violation) arc defined by 
means of the corresponding three-point function 

(p'\JZ m \p)=u(p')T v (q 2 )u(p) , (1) 

where J" m denotes the electromagnetic current and 

r" = 7^1 (<Z 2 ) - 2^<™ V* (« 2 ) + 1 W - 2m ^) ^ Fa (« 2 ) - ^ Vq ^ F * («*) ' W 

with = (p' —p)ix- Here, F\ and F 2 denote the P-, CP-conserving Dirac and Pauli form factors, and m 
is the mass of the nucleon. The last two form factors Fa and F 3 stem from P- and CP-violating terms, 
respectively: Fa denotes the anapole form factor and F3 the electric dipole form factor. In what follows, we 
will only consider the dipole form factor F 3 . The electric dipole moment of the neutron/proton is defined as 
the electric dipole form factor at q 2 = 



A - F 3,n,p(0) ,, 

a n , p - 2m . [6) 



Expanding the form factor in the squared momentum transfer allows one to define an electric dipole radius 
,dF 3 (q 2 



(r 2 ed ) = 6- 



dq 2 



(4) 

<j 2 =0 



Note that similar to the case of the neutron electric form factor, we do not include the normalization of the 
form factor at q 2 = in this definition. 

3. Consider three- flavor QCD in the presence of strong CP violation, parameterized by the constant #0, 
Cqcd = -\Gl v G a ^ + q {1$ -M)q+ Oo^G;„G a ^ (a = 1, . . . , 8) , (5) 

with the gluon field strength tensor G^ v and its dual G a — e^u\aG a ' X(y , q collects the various quarks, and 
is the gauge-covariant derivative. The last term leads to the U(\)a anomaly and is responsible for the 
non- vanishing mass of the rj in the chiral limit. We want to analyze the effects of strong CP violation in 
the appropriate effective field theory, which is chiral perturbation theory. To this end, we treat the vacuum 
angle #0 as an external field and use the appropriate effective Lagrangian for the U(3)l x U(3)r theory. The 
original, systematic construction of the meson Lagrangian for this symmetry can be found in [13]. However, 
for our purposes it is more convenient to adapt the notation given in [7], which is particularly suited for 
the calculation of the electric dipole form factor of the neutron. This formulation itself is partially based 
on [14,15]. As a basic building block, one introduces the external field 6{x) that transforms as 0(x) — > 
6(x) — 2Nfa under isosinglet axial rotations, with Nf the number of active quark flavors. The mesons (the 
eight Goldstone bosons and the singlet field, the 770) are incorporated in a 3 x 3 matrix-valued field U 
that transforms as U — > LUFfi under U(3)l x U(3)r. Since the phase of the determinant of U transforms 
according to lndet U — > lndet U + 2iNfa , one introduces the invariant combination 

6 = 6 -ihidetU , (6) 

which is more convenient for the construction of the effective Lagrangian. The most general effective meson 
Lagrangian to second chiral order, complying with U(3)l x U(3)r symmetry, then reads 

C = -V + V x tr [V AI L> t V^L>] + V 2 tr [x f U + x^] + iV 3 tr [x f U - X U^] 
+ V 4 tr [UV^] tr [U 1 V^?/] + V 5 tr [V^V^tf] , (7) 

with x = 2So (s + ip) and V M f7 = d^U — ir^U + iUl^, s, p, 1^, r M are the standard external sources [13]. 
The Vi are functions of 9. To further analyze these, one makes use of the large- iV c approximation to QCD. 
The use of an expansion in powers of 1/N C in addition to the usual expansion in powers of small momenta 
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and quark masses implies an extension of the power counting scheme. It is convenient to choose the unified 
counting rules according to [14,16] 

p = O (5) , m q = O (<5 2 ) , 1/N C = O (S 2 ) , (8) 

with the pertinent small parameter 5, and m q denotes any of the light quark masses. This allows to expand 
the Vi in powers of 9, with an odd function V3 in 9 and all others even; the expansion coefficient of Vi 
of order n will be denoted by Vj n \ To utilize the effective Lagrangian equation (7), one must fix the 
vacuum expectation value U of U by solving the classical equations of motion. U a can be chosen diagonal, 
expressed in terms of the so-called quark angles. This allows us to write U = \/UqU\/Uq, and choose the 
parameterization 

u = exp {\fli m+ j-/)- (9) 

The coupling constant F for the 770-singlet is in principle different from its octet counterpart due to the 
fact that the unbroken subgroup U(3)v does not exhibit an irreducible representation of dimension nine. 
The resulting effective Lagrangian reads 

£0 = -V + Kitr [W^WU] + (V 2 + BV 3 ) tr [ X (U + £/+)] - iAV 2 tr [U - 

+ i(V 3 - BV 2 ) tr [ X (U - rf)] + AV 3 tv [U + rf] + V 4 tr [C/V ^} tr [uW^U] , (10) 

where the hermitian matrix X nas absorbed further factors of Uo, and A, B are complicated functions of the 
Vi, see e.g. [7]. In leading approximation, they are given by 

Note that 9q = 0(S 2 ), see below. After the vacuum alignment and with our particular choice of U, the Vi 
are now functions of the combination 9q + v^^o/^b- Furthermore, the correct normalization for the terms 
quadratic in the singlet field is obtained by considering all contributions to its kinetic energy from the Vi- 
and V4-terms and demanding the resulting coefficient to equal 1/2. Finally, we can express 9 in terms of 

measurable quantities. For m Ui d <C m s , to lowest order in the quark angles, and neglecting numerically small 

(2) 

corrections to the value of Vq ' [17] this relation reads (in what follows, we set F$ = F n ) 

F 2 M 2 
8V (2) 

In a similar manner, one can construct the most general effective Lagrangian including also the baryon 
octet B up-to-and-including terms of second order in the derivative expansion (we only display the terms 
relevant to our calculation; for details, we refer to [7]): 

£ 0B = itr[B-f[D ll , B]] -mtv[BB] - £^ti[B^ l5 [ Ufll B] ± ] - ^ tv[Bj^ 5 B]tT[u^] 
+ b D/F tr [B [ X+ - iA(U -U^),B] ± ]+b tv[BB] tr [ x+ - iA{U - £/*)] 



+ AAw' w — m tv[BB] +i(w[ 3/1 J + w 13/u —rio)tr[Ba^^F+,B} 



+ w 16/17 tv[Ba^[F+,B] ± ] , (13) 

utilizing U = u 2 , T M = W (<9 M — ir^)u + u(d^ — iZ^) vf] /2, and [D^,B] = d^B + [r M ,S]. Furthermore, 
u M = i [it* (8^ — ir^) u — u (<9 M — il^) ut] , x+ = ux^u + u^x u ^ ancl F+ v incorporates the field strength tensor. 
Also, rh is the octet mass in the chiral limit, F, D are the conventional axial coupling constants, A is an 
isosinglet axial coupling, and the b / D / F are the low-energy constants (LECs) related to the leading-order 
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Fig. 1. Feynman diagrams contributing to the electric dipole form factor of the nucleon. Diagrams (a), (b) contribute at tree 
level, (c)-(h) are of the second and (i)-(l) of the third chiral order. Filled circles are second-order mesonic vertices, squares 
and diamonds represent vertices generated by the first- and second-order meson-baryon Lagrangian, respectively. CP-violating 
vertices are denoted by ®. 

symmetry breaking. The LECs Wi parameterize the coupling of the singlet field to the baryons at second 
order (we use the abbreviation w[ — u>io + 3u>i2/2 compared to [7]), except for W16/17 which are the 
conventional magnetic moment couplings. 

4. We now have assembled all pieces to calculate the dipole form factor of the neutron (and the proton) 
to third order in the chiral expansion. In Fig. 1 we show the corresponding Feynman graphs. These are 
tree graphs (a), (b), one- loop graphs starting at 0(5 2 ) (c)-(h), and one- loop graphs with one insertion 
from the dimension-two chiral meson-baryon Lagrangian starting at 0(S 3 ) (i)-(l). To obtain a consistent 
power counting, the loop diagrams are evaluated using the so-called infrared regularization scheme [18]. In 
this covariant framework, loop graphs resum an infinite string of higher-order corrections, associated with 
the 1/m insertions in the heavy-baryon propagator. In the following, we will expand these contributions 
up-to-and-including 0(S 3 ). It can be shown that the graphs (g)-(l) do not contribute at this order or cancel 
each other exactly. For later purpose, it is convenient to decompose the form factor into tree level and loop 
contributions, F3 — F^ rcc + F^ 009 . After renormalization, the tree level part for the neutron is given by 



T-itrcc 



-16me0Q 



-w 
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48 



T/(1)t/(2),„ 



(14) 



where w[ 3 is the finite part of w' 13 . The third-order loop result for the neutron electric dipole form factor 
follows as 



(D + F)(b D 




1-ln 



o-r In 



Ml 
1 



ov In ■ 



1 ir q 2 - 2Ml 



2 m\/ ' —q 2 



a K -l 



it q 
2~ 



0V - 1 

2M% + 4mAm S jv 



arctan 



arctan 




where, in the space-like region, q 



-Q 2 < 0, Am SJV 



m, which is of order i5 2 , and u^/ K 



1 — AM 2 , K / q 2 . The second-order result is given by omitting the last term of the pion and kaon part in the 



expression above. In the limit q 2 — > 0, one can readily deduce the expression for the electric dipole moment 
of the neutron, 



d n = dT° + <4°° P , dT c = -8 e 9 Q< 3 + F 2pt m V^VJVa 



± TT^VO 



d l ° op 



eVr 



(2) 



T 2 Fj 



(D + F)(b D +b F ) 1 + ln^ 



2to 



-(D-F)(b D -b F ) ( 1 + ln ^ 



M£ ttMk tt Am E jv 



2m 



+ 



M 



K 



(16) 



The expression for the tree result was already derived in [7] , however we obtain an additional factor of three 
for diagram (b). Note that we have kept the dependence on the scale of dimensional regularization /i in 
the loop expression in order to give an error estimate by varying the scale when performing the numerical 
evaluation. Dropping all terms besides the chiral logarithms reproduces the result of [9,7]. The next-to-leading 
order correction in the pion-loop contribution was also obtained in [19] in a relativistic loop calculation. We 
have also calculated the corresponding form factor of the proton to achieve an isoscalar/isovector separation 
in the two- flavor case as discussed below. We refrain from giving the corresponding formulas here [20] . From 
the neutron and proton form factors one can derive analytic results for the corresponding radii, the one for 
the neutron reads 



( 



rl d ) r 



2eV r 



(2) 



r 2 Fi 



— (D — F)(pD — bp) 



Ml 



5tt 
4M K 



7rm Am 



SAT 



2M\ 



(17) 



5. Before showing results, we must fix parameters. We use D — 0.804 and F — 0.463 from hyperon /3- 
decays [21], bp — 0.066 GeV -1 and bp = — 0.209 GeV -1 from the leading-order analysis of the baryon mass 
splittings, F ff = 92.4 MeV as well as V {2) = -5 • 10~ 4 GeV 4 and V^ 1] = 3.5 • 10~ 4 GcV 2 from an analysis of 
r] — i]' mixing [17]. The same analysis also shows that the correction to F^/F^ = 1 is of order 1/N C , so that 
we use F = F^ based on the counting rule in Eq. (8). 

First, we consider the neutron dipole moment. Following the arguments of [7] based on large- N c counting, 
the ratio of diagrams (a) and (b) is given by 



,tree(b) 



,trcc(a) 
d n 



144 



F 4 M 2 



0.38 . (18) 



This ratio will serve for an error estimate of the tree level contribution. Furthermore, the Wi6-term of the 
meson-baryon Lagrangian, describing a magnetic coupling to the baryon fields, is similar to the structure 
proportional to w' 13 / w\ 3 except for the additional 75. However, the coefficient wiq is of O(N^), whereas 
w\ 3 is of (D(N®), which suggests the assumption |w[ 3 | < \wie\. Utilizing the third-order calculation of the 
baryon electromagnetic form factors [22], we find wis = 0.40 GeV -1 . Assuming further |w[ 3 | = |wi6|/3, we 
obtain a bound on the dipole moment, 



|(4 roe | ~ (2.9 ±1.1) x 10 -iD 6> o ecm . (19) 

The numerical evaluation of the loop corrections to the electric dipole moment of the neutron is much more 
straightforward since it does not involve unknown parameters. We find 

d l ° op = S.Otl l x 10 -16 O e cm , (20) 

varying the renormalization scale fi between M p and m~ (note the much larger scale variation in [19] that 
induces an even more pronounced uncertainty). However, it was shown in [23] that the corrections to b D , 
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Fig. 2. Loop contribution to the electric dipole form factor of the neutron as a function of Q 



-q to third (red, solid 



line) and second (blue, dot-dashed line) order, using /i = 1 GeV. The boundaries of the error bands due to variation of the 
rcnormalization scale /i between M p and m= are represented by red dashed or blue dotted lines, respectively. 

&f due to fourth-order effects are large, so that the central value for the loop contributions changes to 
d\°° p = —5.1 x 1CP 16 #o e cm. Taking this latter value, the loop corrections slightly dominate and we can 
deduce a lower bound for the theoretical estimate of the electric dipole moment, 



jthco I 



> 1.1 x 10~ 16 6> o ecm 



(21) 



Together with the experimental upper bound on the electric dipole moment of the neutron, this finally yields 
|0 O | < 2.5 x 1(T 10 . (22) 

We remark, however, that this bound is very sensitive to cancellations between the tree and loop contributions 
and based on large- N c arguments, and thus should be considered with caution. 

Next, we discuss the loop contribution to the neutron dipole form factor. It is shown in Fig. 2 in units of e 
times the dimensionless parameter 9q . The full result including third-order corrections is represented by the 
solid line, whereas the dashed line contains second-order contributions only. The error bands of both lines 
are obtained by varying the renormalization scale between M p and ms. again. The renormalization scale for 
the central curves has been set to the neutron mass /i = to. As can be seen from this figure, the relative 
size of the third-order corrections compared to the leading-order terms is small for reasonable values of Q 2 . 
In fact, the change due to the variation of the renormalization scale clearly exceeds the change resulting 
from third-order corrections. Nevertheless, the slope of the form factor for small values of Q 2 is altered 
significantly by the next-to-leading-order corrections, as it is also borne out by evaluating Eq. (17). We find 



-20.4 [1 - 0.67 + 0(S 2 )] (9 efm 2 



(23) 



and a similar effect is found for the proton. This large correction can be traced back to the additional factors 
of 7r in the pion loop contribution in Eq. (17). Such factors also appear in the analysis of the radii of the 
isospin- violating nucleon form factors, see [24]. We also note that the contribution from kaon loops is much 
smaller than the one from the pions and therefore, the radius is almost entirely of isovector nature. 

To compare results from two-flavor lattice QCD at unphysical quark masses with predictions from chiral 
perturbation theory, it is necessary to perform an extrapolation of the analytic results in the pion mass. To 
this end, we have to express Ato^jv in terms of M n , Mk and account for the fact that Mr- itself depends 
on M w . To lowest order this dependence is given by 



M 2 K = M 2 K 



1 



-Mi 



(24) 
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Fig. 3. Plot of d n p as a function of the pion mass to third (red, solid line) and second (blue, dot-dashed line) order, 
with fi = 1 GeV. The boundaries of the error bands due to variation of the renormalization scale fi between M p and m= are 
represented by red dashed and blue dotted lines, respectively. The data from two-flavor lattice QCD [4] are shown by the black 
squares. 

where Mk denotes the mass of the kaon in the SU(2) chiral limit. Numerically, one has Mk — 0.484 GeV. 
Special care has to be taken furthermore to relate 9o to 9o, as Eq. (12) has been derived assuming m u ,d <§C m s , 

(2) 

as well as a certain hierarchy of terms relying on the numerical value of Vq . Giving up those assumptions 
changes the relation between 9q and 9q to 



In Fig. 3 we show the resulting pion (quark) mass dependence of the loop contribution to the electric dipole 
moment in comparison to the available data points from two- flavor lattice QCD [4]. It is interesting to see 
that the third-order calculation reproduces the trend of the lattice data (the order of magnitude and the 
global sign) even without the unknown tree contribution. However, only below pion masses of the order of 
500 MeV the third-order corrections are sufficiently small for a stable chiral extrapolation. 

6. Finally we consider the matching of our three-flavor results to the two-flavor representation. This not 
only allows to extend the findings of [10] to higher order, but also gives an explicit representation for some of 
the LECs appearing in that paper. Such matching relations between three- and two-flavor versions of chiral 
perturbation theory have been studied extensively in the meson sector (see the most recent comprehensive 
results [25,26] and references therein), and also first results for the meson-baryon sector exist [27,28]. 

The leading-order SU(2) formulas for the electric dipole moment of the nucleon have been worked out 
in [10], within the framework of heavy-baryon chiral perturbation theory. The authors employ a decompo- 
sition of the electric dipole form factor of the nucleon according to 




(25) 




(26) 



where F^°\q 2 ) is the isoscalar and Fg 1 ^^ 2 ) is the isovector part of the electric dipole form factor. The 
relations between the corresponding isoscalar and isovector electric dipole moments d /i and those of proton 
and neutron are given by 
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By matching our SU(3) result to the two-flavor theory, we can extend the expressions given in [10] by one 
order in the chiral expansion, 



d = d Q + 



8n 2 F v 



3tt M„ 
4 m 



/ , Ml SirMv 

d\=di+ — 1 + In - 

8n 2 F v V fi 2 4 m 



(28) 



where the d /i parameterize the short-distance physics that now includes kaon, i]s, and 770 loops. Explicitly, 
the matching relations for these two constants read 



<V = \ (dp™ + 4°° P ) ± J (^ r °° + ^° P ) + ©(M^, M^J , 



d$ ee 



<? p oop 



48Vy'V r 



(1)t/(2) 







2 +12^ol^^ 2 



-Wl4 



Off r/ 



jloo 



(2) 







2 TT^ 
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c (£» + 3F)(6r, + 3M U + m 
6 \ [i 



Jtree o 5 f W 13 48V 3 (1 V (2) 



2 



(D - F)(6 D - b F ) + 



(D - 3F)(b D - 3b F )\TrM K 



3^3 



m 



F 2 



-(2D + 3 A) (2b D + 3bo + 6V3w' w ) 



c?° op = 



(2) 



(D-F)(6d-H(1 + ^ o 

A 4 



2m 



(29) 



where we have used the SU (2) chiral limit relation M m = 2Mk/V3- The new terms in the SU(2) expressions 
Eq. (28) are the ones linear in the pion mass. We note that the corrections ~ M n are again sizeable. 
Furthermore, it was argued in [10] based on dimensional analysis that the CP-violating pion-nucleon coupling 
g-KNN has to be of order 



g-xNN — O 



OqM 2 
mFif 



(30) 



The previous results allow us to obtain a matching relation for q^nn • Utilizing the standard matching 
relation for the two- and three- flavor axial couplings, gA = D + F + 0(<5 2 ), we obtain for g^NN 



g,NN^^-(b D + b F )+0(S 4 ) , 



(31) 



which numerically yields g^NN — — 0.03[— 0.05] 9 , where the number in the brackets refers to the values of 
b D j F from [23]. This confirms the dimensional estimate equation (30) and reproduces the result originally 
obtained in [8]. However, our representation in terms of symmetry-breaking LECs only is certainly more 
compact and can be systematically improved by going to yet higher orders in the chiral expansion. 

7. In this Letter, we have analyzed the neutron (and proton) electric dipole form factor and its moment in 
the framework of covariant U(3)l x U(3)r baryon chiral perturbation theory, extending earlier calculations 
by one order in the power counting. The latter includes the expansion of QCD in the number of colors. 
While the third-order corrections to the neutron electric dipole moment are small at the physical pion mass, 
their contribution makes the trend of the quark-mass expansion agree with the existing data from full lattice 
QCD, although those still employ pion masses too high for a safe chiral extrapolation. More lattice data 
at smaller pion masses will allow one to analyze strong CP violation as encoded in the neutron electric 
dipole moment in more detail. From comparison with experimental limits, we have given an upper bound 
on the vacuum angle. The electric dipole radius, given entirely by chiral loops, receives a large correction at 
next-to-leading order. Finally, we have also given matching relations for our three-flavor representation to 
the SU(2) case. 
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